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Summary. We provide a variational description of any Liouville (i.e. volume 
preserving) autonomous vector fields on a smooth manifold. This is obtained via 
a "maximal degree" variational principle; critical sections for this are integral 
manifolds for the Liouville vector field. We work in coordinates and provide 
explicit formulae. 

Introduction 

It is well known that Hamiltonian dynamics preserves the volume in phase space 
(Liouville theorem); they are thus a prominent example of incompressible (or 
volume-preserving) dynamics. 

However, other dynamics which preserve the volume, or more generally a 
measure, in phase space without being necessarily Hamiltonian, are also of 
physical interest. These are called Liouville vector fields (see below for a pre- 
cise definition). The motion of an incompressible fluid is of course described by 
volume-preserving dynamics, and is a prominent physical example of Liouville 
dynamics. 

Among relevant classes of volume-preserving dynamical systems which are 
not necessarily hamiltonian, we mention Nambu mechanics [15] and reversible 
dynamics [10]. 

It has long been known that several subclasses of Liouville dynamics share 
important features with Hamiltonian one; for example, statistical mechanics can 
be properly based on Nambu dynamics (this was actually the motivation for it 
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[15]), and reversible vector fields share many results of the perturbation theory 
of Hamiltonian systems, including KAM theory [10]. 

It is thus entirely natural, and justified by physical relevance, to wonder what 
is the extent to which Liouville dynamics shares the structures and properties 
of Hamiltonian one. 

The question has been studied by several groups; we mention here in partic- 
ular the work of [11], extending to Liouville dynamics several geometric struc- 
tures - including the Poisson bracket formulation - of Hamiltonian dynamics, 
and clarifying several points in the geometrical formulation of Liouville dynam- 
ics; see [7, 12] for further developements in this direction (also connected with 
Nambu-type systems). 

Hamiltonian dynamics can be characterized in terms of a variational prin- 
ciple, and this is instrumental to many results in hamiltonian theory; it is thus 
natural to investigate if some kind of variational formulation is also possible for 
general Liouville vector fields (in general, a standard variational formulation is 
not possible). 

The purpose of this note is to answer this question in the positive: we find 
that any Liouville vector field on the phase space manifold P can be described as 
the unique properly normalized characteristic vector field of a maximal degree 
variational principle (these notions will be defined below) on the extended phase 
space M = R x P; this applies in particular to autonomous Liouville vector field 
on P, in which case the R factor can be thought as the physical time. 

It should be stressed that, although in this note we work in local coordi- 
nates, an intrinsic discussion would be possible and is actually given elsewhere, 
leading to more general results [5]. Such an approach, however, requires to use 
the Cartan theory of exterior differential systems; this theory is well known in 
differential geometry and in the geometric theory of differential equations, but 
we believe the results given here can be of interest to all physicists working 
on incompressible dynamics, which in general are not necessarily familiar with 
Cartan theory. 

Luckily, the case of Liouville fields can be dealt with at the much simpler level 
considered here, i.e. working in coordinates and involving only basic notions of 
differential geometry and rather elementary mathematical analysis. 

Acknowledgement. The work of GG has been supported by "Fondazione 
CARIPLO per la ricerca scientifica" under the project "Teoria delle pertur- 
bazioni per sistemi con simmetria" . 

1 Liouville dynamics 

1.1 Liouville vector fields on phase space 

Let P be a smooth and oricntablc A^-dimensional manifold, from now on called 
the "phase space"; we will, as customary, denote by A(P) the set of differential 
forms on P, and A k (P) C A(P) denotes the set of forms of degree k. Recall that 
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all volume forms on P (also referred to as measures) are equivalent; we choose 
a reference volume form £1 on P. 

Let X be a vector field on P; we say that X is a Liouville vector field (with 
respect to the measure £1) if it preserves £1, i.e. if Cx(Q) — 0, with C the Lie 
derivative. 

As C x {a) := {X _lda) + d(X _ia) and d£l = 0, X is Liouville w.r.t. £1 if 
and only if X _l £1 is closed; locally this means that there exists a (N-2)-form 7 
such that 

A_l£l = d 7 . (1) 

Thus X is Liouville with respect to £1 iff for any neighbourhood iCP, eq.(l) 
is verified for some 7 G A JV ~ 2 ( J 4), i.e. iff A_l£l is exact on A. 

We say that A is globally Liouville with respect to £1 iff (1) is verified for 
some 7 G A Ar_2 (P), i.e. iff A_l£l is exact on P. 

In the following, for ease of discussion and notation, we will deal with the 
case where X is globally Liouville; however it is clear that for the sake of local 
considerations the two cases are equivalent. 

Note that if X and £1 are given, 7 is not uniquely defined by (1) (we can 
always add a closed form 71); on the other hand, if £1 and 7 are given, then (1) 
uniquely defines X. 

We have chosen a given volume form £1; for any different choice £1 of the 
volume form we have £1 = pCl, with p a nowhere vanishing (positive, to preserve 
orientation) function p G A°(P); moreover if F_l£l = d7 and £1 = pil, then 
X := pY satisfies (1) with the same 7. Hence whenever we have a vector field 
Y which is Liouville with respect to a measure £1, we can reduce to considering 
X and an exact volume form. Thus, from now on we will assume £1 = da, for 
ease of computation. 

1.2 Enlarged phase space 

It is convenient to consider the "enlarged phase space" M :— R x P, where the 
R factor corresponds to the time coordinate t; we denote the reference volume 
form in M corresponding to £1 G A N (P) by £1 M := dt A fl G A N+1 (M). Then 
the dynamics defined by X in P corresponds to the flow of Z in M, where 

Z := d t + X ; (2) 

this satisfies, by construction, Z _ldt = 1. 

We consider now a form G A N ^ 1 (M) built from the a defined above and 
from the form 7 associated to X via (1); this is defined as 

■d := a + dt A 7 ; (3) 

it is immediate to check that dd is nowhere zero. 

We now prove that the form •& defined in (3) defines a unique vector field Z 
on M via 

Z_ldi9 = ; Z_idt = 1 . (4) 
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Note that Z _ldt ^ implies that the vector field Z has a nonzero component 
along d t . 

To see that (4) defines a unique field, it suffices to note that d§ is a TV-form 
in the manifold M of dimension N + 1: in this case its annihilator Af(d{}) (the 
set of vector fields Y on M such that Y _l di? = 0) is a one dimensional module 
over A°(M). Obviously the first of (4) amounts to the requirement Z G jV(d#). 

Remark 1. More precisely, any form a G A JV (M) can be written in local 
coordinates (x°, x 1 , x N ) (we use the notation <9 M = djdx^ and £ = x°) as 
A M (x)(9 /i —I^m); it is easy to check that the vector fields in M{a) are then 
written as Y = with f»{x) = F(x)A»(x), where F G A°(M). 

Note next that Y G A/"(dt?) (and F ^ 0) is necessarily such that Y _\dt ^ 0, 
just by the form of The second equation in (4) is therefore just a normalization 
condition, selecting a unique vector field Z out of the one dimensional module 
Af(dtf). 

Actually, the vector field Z defined by (4) is just the vector field Z defined 
above in (2). Given that Z is unique, to prove that Z = Z it suffices to check 
that Z satisfies (4). 

It is obvious that Z _\dt — 1, as already remarked. As for the first of (4), 
we have, using da = ft and (1), 

Z_\d9 = X_M + (-l) N (Z_\dj) Adt + (-i)( 2JV - 1 )d 7 
= +d 7 + (-l) JV (X_ld7) Adi - d 7 
= (-lf(Xjd7)Adf = (-lf[U(Xj!])] = 0. 

This completes the proof; we summarize our results as follows. 

Lemma 1. The equations (4), with ■d given by (3), select uniquely the vector 
field Z given in (2), where X satisfies (1). 

2 Variational setting 

We want now to prove that Z , and hence A, can be given a variational char- 
acterization. This will be based on a suitable fibration tt : P — > Q of the 
phase space over a smooth if-dimensional manifold Q, which induces a fibra- 
tion it : M — > B of the extended phase space, ir := id x n. In order to do this 
we will have to introduce the concept of maximal degree variational principle, 
and of the associated characteristic field. 

2.1 Variational principles 

Consider a general smooth bundle tt : E — > B, with E a smooth manifold of 
dimension n and B a smooth manifold of dimension k, with 1 < k < n (in the 
next section we will specialize to the case E = M = H x P). We assume the 
fibers 7r _1 (x) arc parallelizable. 
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We denote by r(7r) the space of smooth sections ip : B — > E of this bundle, 
and by V(ir) the space of vector fields on E which are vertical for 7r, i.e. tangent 
to fibers n^ 1 (x). Given a form a G A(£ l ), we denote its pullback by a section 
<p S r(7r) as <p*(a). 

Let D C 5 be a domain (i.e. a closed compact manifold with boundary) in 
B, and r] G A k (E); we define a functional I : T(tt) — > R by 

Ifc) := / ^(r?) . (5) 

A vector field V G V(7r) obviously induces an action on r(7r); this results in 
turn into an action on I (if). More precisely, denote by ip 3 the local flow of V 
on E, and consider a section ip G r(7r); the flow in r(7r) originating from <p is 
the one-parameter family of local sections f s = tp s (fo) : = "0s <Po- We define 
the variation of I under V as 



(6 v I D )(f) := ^ 



D 



(6) 



s=0 



We denote by Vd(7t) C V(tt) the set of vector fields vertical for it which 
vanish on the fibers over 3D. We say that ip G r(7r) is extremal (or also critical) 
for J defined by (5) if and only if 8yl{ip) — for all V G Vd(7t). 

It is well known that, equivalently, ip is critical for / defined by (5) if and 
only if 

ip* (V _l dr]) = W e Vd(tt) . (7) 

Given the complete equivalence of this condition with the previous definition 
[9, 16], we can take (7) as the definition of critical section. 

It will be convenient, as our considerations will be local, to use local coor- 
dinates in E adapted to the fibration it : E — > B; that is, we use coordinates 
(x 1 , ...,x k ;y 1 , y p ) where k + p = n, (x 1 , ...,x k ) will be coordinates in B, and 
(y 1 ,...,y p ) are vertical coordinates spanning the fibers tt^ 1 (x). We write, for 
ease of notation, 

d t :=d/dx l , 8 a := 8/dy a ; i=l,...,k, a=l,...,p. (8) 

The image of a section ip e r(7r) is then locally the graph of a function u : R fe — > 
R p , i.e. f — {x, y = w(x)}. In this coordinates, ip*(dy a ) — (du a / dx l )dx l . 

The determination of equations describing critical sections for I is very well 
known, and we just recall it to fix notation. Any V G V(tv) is written in these 
local coordinates as V = f a d a ; the set Vd(tt) is identified by the condition that 
f a = on tt-^OD). Thus (7) requires f*[f a {d a _ldi?)] = for arbitrary f a 
(vanishing on n^ 1 (dD)); hence wc need that p*(d a _ldi9) = for all a — 1, ..,p. 

It is natural to introduce the fc-forms ^ a '■— da—ldfl; the section ip is critical 
for / if and only if 

f*(V a ) - a = l,...,p. (9) 

Remark 2. It should be noted that the fibration ir : E — > B is nontrivial 
for 1 < k < n — 1; however, the case k = n — 1 will provide trivial equations 
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(as discussed in [5]); thus we have to require 1 < k < n — 2. The case k = 1 
corresponds to the well known case of one independent variable (as in standard 
Hamilton dynamics) , and it is well known that in this case the variational princi- 
ple identifies a vector field on E, i.e. yields a system of ODEs for ip. For k > 1, 
the variational principle yields a system of PDEs for ip; in local coordinates 
as those introduced above, this is a system of p equations for the p functions 
^"(x 1 , ...,x k ) depending on k independent variables. 

Remark 3. Note that the fibration tt : E — > B is to a large extent arbitrary; 
basically it amounts to a choice of which variables should be considered as inde- 
pendent ones and which as dependent in formulating the variational principle, 
and has no intrinsic meaning. Later on, when we specialize to£ = M = RxP 
and we study a dynamics (see sect. 3), we will require that the time is kept as 
one of the independent variables: this means that B = RxQ, and the fibering 
7r : M — ► B is induced by a fibering n : P — ► Q as anticipated above. Note also 
that the variational principle is always formulated over a domain D C B, i.e. in 
local terms. 

2.2 The maximal degree case 

In view of remark 2 above, we say that for k = n — 2 we have a maximal degree 
variational principle. We are interested in this case: we will show that these 
variational principles still identify a vector field on E, i.e. a system of ODEs 
associated to it. 

As k = n — 2, we have p — 2; we will write z = y , w = y 2 to avoid a 
plethora of indices. We write u> — dx 1 A ... A dx k for the reference volume form 
in B\ the reference volume form in E will of course be 7r*(w) A dz A dw; in the 
following we will write, with a slight abuse of notation, oj for w*(u>). In this case, 
when studying the variational principle defined by (5), (6), we should consider 
(3 := dn e A n ^ 1 (E); we can always write any (3 £ A n ^ 1 (E) in the form 



with fi = 1,2, A^,f,g smooth functions of (x, z,w), and := _l uo. 

In the following, we will assume that the vector A = (A 1 , A k ) is not 
identically zero; in this case we say that the maximal degree variational principle 
defined by n is proper. 



In the present notation, we choose d z and d w as generators of V(n), i.e. 
*i = d z -i/3, * 2 = du,— 1/3. With ip e T(ir), we have 



/3 = EjUi A? [u w Adz A dw] + 

+ (-l) k f [uj A dw] + (-l) k+1 g[uj Adz] 



(10) 



*i = (-l) fe_1 [A"(w w Adw) + (-l) k gu;] ; 

*2 = (-l) fc [4"(w w Adz) + {-l) k .fuj} . 

<p*(^i) = <fi* [A^(dw/dx») - g]u; 

<P*(y 2 ) = -if* [A^{dz/dx^) - f]u . 



(11) 
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Requiring the vanishing of both ip*{^> j) for j = 1,2 means looking for solu- 
tions of two quasilinear first order PDEs, i.e. 

p*[C Y {z) - f] = ; <p*[C Y (w) - .9] - (12) 

(with obviously Y = A^d^, and Ly the Lie derivative). Note that as A ^ in 

(10) , we are guaranteed 7^0. 

The relevant property is that the equations can be written in terms of the 
action of the same (nonzero) vector field Y , or more precisely [1] in terms of the 
(non vertical, as Y ^ 0) vector field W = Y + fd z + gd w on E, i.e. 

W = E^(x;z, W )^r + /(*;*,«0^ + 9tez,v>)£;- (13) 

Indeed, as well known (see e.g. chapter II. 7. G of [1]), the Revalued function 
u(x,t) = (z(x,t),w(x,t)) is a solution to the system of quasilinear PDEs (12) 
if and only if its graph is an integral manifold for the associated characteristic 
system 

dx> i /ds = A> 1 , dz/ds = f , dw/ds^g , (14) 

i.e. for the W given above (the characteristic system (14) is often written in the 
so-called symmetric form, i.e. as dx 1 /A 1 = ... = dx k /A k = dz/f = dw/g). 

We recall that if X is a vector field on E and the submanifold S C E is 
such that X(x) G T X S for all x € S, we say that S is an integral (or invariant) 
manifold for X; if S is invariant for X and one-dimensional, we also say it is an 
integral curve of X. 

It is thus entirely natural to call the W given by (13) the characteristic 
vector field for the maximal degree variational principle on tt : E — > B defined 
by a form n G A n ~ 2 (E) such that dr] — (3 is given by (10). 

Note that the above computations imply that for given 77 we have a one 
dimensional module of characteristic vector fields (all differing by multiplication 
by a nowhere zero smooth function), defining a unique direction field on E (see 

[11) - 

Summarizing, with the above discussion we have proved that: 

Theorem 1. Let tt : E — > B be afibration of the smooth n-dimensional manifold 
E over the (n — 2) -dimensional manifold B, and let r) G A n ~ 2 (E) be such that 
dr) = (3 is written as in (10). Then the section ip G T(7r) is critical for the 
maximal degree proper variational principle defined via the functional (5) if and 
only if it is an invariant manifold of the characteristic vector fields for the 
variational principle. 

Corollary. Given an integral curve Y of a characteristic vector field for the 
variational principle defined by (5) and a critical section ip of it, either T Dip = 
or T C tp. 

The corollary implies that we can describe the /c-dimensional critical sections 
<p for a variational principle as the union of integral curves of the characteristic 
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vector field W for the same variational principle passing through a suitable 
submanifold po C <p of dimension k — 1 . 

Note that if one of the A^, say A 1 , is nowhere zero, we can divide this out 
from W, and obtain a vector field of the form 

Z = di + X (15) 

with X-ldx 1 = 0. 

We would also like to remark that W is tangent to sections ip s r(7r) such 
that if*(^i) = = if* (^2), see above and [1]; on the other hand, <£*(*») = 
means that \&i vanish on vector fields tangent to ip, hence vanish if evaluated 
on W, i.e. W_l*j = 0. We have thus proven that: 

Lemma 2. T/ie W identified by (13) satisfies W_l*i = = W_l* 2 - 

3 Variational principle for Liouville dynamics 

In this section we apply the discussion on maximal degree variational principles 
developed in the previous section, to the case of Liouville dynamics. 

We consider, with the notation introduced in section 1, £ = M = R x P, 
and B = R x Q; note that n = N+ 1 and k — K + 1. We stress that the fibering 
7r : M — > B is naturally induced by a fibering n : P — > Q, i.e. with obvious 
notation n = idxn (see also remark 3 above). 

Accordingly, we will use local coordinates adapted to the double fibration of 
M; we call these (t; (; k ; y 1 , y 2 ) , where t is the coordinate along R, (£*, ...£ fe ) 

are coordinates in Q - hence (t; ^) are coordinates in M - and y a are "vertical" 
coordinates on the fibers 7r _ 1 (t ; £) = 7?^" 1 (^) . The setting of the previous section 
is recovered via x 1 :— t, := 

In order to recover the Liouville vector field X satisfying (1) as the vector 
field identified by a maximal degree variational principle, it suffices to choose as 
T) the form # defined in (3). Note that in this case, using the notation (10), we 
have indeed that A 1 is nowhere zero. 

With local coordinates (x, z, w) as in the previous section, write V\ = d z , 
V2 = d w and consider a vector field X which is nonzero and non vertical; hence 
Vj _l (X _l dn) — for j — 1,2 means that \ X _ I dr] does not contain dz or 
dw factors. However, this is impossible unless X _idi] — 0. Indeed, \ S A fe (M), 
hence \ — ^( x i z ; w)dx x A ... A dx k \ this cannot be obtained by x = X _i(3 if X 
is independent of <9 Z and d w . In other words, we have proven the 

Lemma 3. Let X be a non-vertical vector field for the fibration tt : M — ► B, 
where dim(M) = dim(S) + 2. Then V _l (X _l dr/) = /or a// V e V(7r) implies 
- and is thus equivalent to - X _idrj = 0. 

Note that this applies to W provided this is non vertical, i.e. provided the 
A^ identifying dr/ - see (10) - are not all identically vanishing. However, as 
remarked above, in the case we are discussing A 1 ^ at all points, and this 
condition is hence satisfied. We have thus proven the following: 
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Theorem 2. The characteristic vector field W for the variational principle 
defined by rj satisfies W —idn — 0. 

Recalling that Af(drj) is one dimensional (see also remark 1 for its explicit 
description) we also have the 

Corollary. Any vector field V such that V —idn = satisfies V = fW for 
some smooth function f : M — > R ; with W the characteristic vector field for the 
variational principle defined by rj. 

It follows from lemma 3 that we can redefine proper maximal degree varia- 
tional principles and their characteristic vector fields of a maximal degree vari- 
ational principle in a coordinate-independent manner as follows. 

Definition 1. The maximal degree variational principle defined by n on the 
fiber bundle ir : M — > B is proper if there are vector fields V\, V2 G V(ir) such 
ttat[Vi_l(V2-ld»7)]^0. 

Definition 2. A vector field W on M satisfying W —idn — is a characteristic 
vector field for the maximal degree proper variational principle on n : M — > B 
defined by n. 

Using lemma 3, the discussion of the previous section and that of section 1, 
when compared, show that the vector field Z defined in (4) is a characteristic 
vector field for the variational principle defined by z9, see (5); note that imposing 
W _l dt = 1 selects among the characteristic vector fields the one with the proper 
normalization, i.e. imposes W = Z. 

This provides therefore a unique characterization of the Liouville vector field 
X in terms of the variational principle on the bundle ir : M — > B defined by the 
form 1?. We summarize our result in the following 

Theorem 3. Let X be a vector field on the phase space manifold P, globally 
Liouville w.r.t. the volume form — da and satisfying (1). Then the vector 
field Z := d t + X on the enlarged phase space M = R x P is a characteristic 
vector field for the variational principle on ir : M — > B defined by the form -Q in 
(3), and is uniquely selected by the normalization condition Z _ldt = 1. 

Obviously, our discussion could as well be read in the opposite direction: 
consider M — R x P where P is TV-dimensional and orientable; then, given any 
form 7 G A N ~ 2 (P), our construction provides a unique vector field X on P, 
which is obviously guaranteed to be Liouville w.r.t. the volume form f2 := da. 

This is similar to the situation met in Hamilton dynamics on a symplcctic 
manifold (S, oj), where the relation between the symplectic form ui G A 2 (S), the 
hamiltonian function H G A°(S) and the hamiltonian vector field X is given by 
X _\lu = dH. On the one hand, any given H generates through this relation 
a unique vector field X. On the other hand, one can assign the vector field X 
and look for functions H (obviously depending on X itself) generating such a 
vector field; as well known this also provides a variational characterization, via 
Hamilton's principle, for the vector field X. 
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Remark 4. Finally, we mention that here we have dealt with vector fields 
which are globally Liouville, see sect. 1.1; however, similar results hold in the 
case of locally Liouville vector fields [11], with obvious modifications: e.g. the 
forms 7 and a are not necessarily globally defined, and one has to work chart by 
chart. Similarly, to avoid chartwise discussion we dealt with the case of globally 
exact volume form, ft = der, but this restriction is inessential, as discussed at 
the end of sect. 1.1. 

4 Examples — I 

In this section we want to consider three general classes of Liouville dynamics, 
i.e. Hamilton, Nambu, and Hypcrhamilton dynamics; specific equations will be 
considered in next section. 

Note that, in view of our results (see in particular theorem 2) the interesting 
object is the form while the choice of the fibration ir : M — > B is inessential 
and has no intrinsic meaning (see remark 3); we will thus provide, in this and 
the next section, just d in order to describe how the considered systems are 
dealt with in our formalism. 

4.1 Hamilton dynamics 

As well known, any hamiltonian vector field is also Liouville. Let us describe how 
this is identified by a maximal degree variational principle (beside the standard 
minimal degree variational principle based on the Poincare-Cartan one- form). 

Let P be a smooth manifold of dimension n = 2m, equipped with a sym- 
plectic form u; we write £ = (l/(m — l)!)(a))( m_1 ) (this is obviously an exterior 
power). Choose a form p e A X (P) such that locally ui = mdp. 



The smooth function H : P — > R defines the hamiltonian vector field X by 
X_iuj = dH. On the other hand fi = (l/m!)(w) m = (l/m)w A (. Thus 



and X satisfies (1) with 7 = H(. 

It follows immediately that the corresponding maximal degree variational 
principle is based on the form $ E A 2m_1 (M) (recall M = R x P) given by 



the Hamilton equations are readily recovered from this. 
4.2 Nambu dynamics 

Nambu dynamics [15] encountered a renewal of interest in recent years, see e.g. 
[8, 17]; see [13] for a discussion of it in terms of forms and Cartan ideals. It is 
well known that Nambu dynamics is also Liouville, and that in general it cannot 
be described in terms of a standard (i.e. degree one) variational principle. 



X_M = (X_iuj)A( = dHA( 



(16) 



7? = pA( + H( Adt = (p + Hdi) AC 5 



(17) 
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An intrinsic definition of Nambu vector fields is as follows: consider a smooth 
n-dimensional manifold P with volume form fl. Then the vector field A on P 
is Nambu if there is a choice of n — 1 smooth functions Hi : P — > R (i — 2, n) 
such that 

dH 2 A...AdH n := x = A_lft . (18) 
We have immediately that \ = with e.g. 

7 = H 2 dH 3 A ... A dH n . (19) 

We also write e.g. a — x x dx 2 A ... A x n , which yields da — 0. (For both 7 and 
a one could use a different permutation of indices). 

With these, 8 is readily recovered, see (3), and hence - for any A - we have 
determined the maximal degree variational principle defining the Nambu vector 
field X. 

4.3 Hyperhamiltonian vector fields 

Another special class of Liouville vector fields is provided by hyperhamiltonian 
vector fields, generalizing Hamilton dynamics and studied in [4, 14]; these are 
based on hyperkahler (rather than symplectic) structures [2] . 

In this case, one considers a riemannian manifold (P, g) of dimension p = AN, 
equipped with three independent symplectic structures u> a (a = 1,2,3); to a 
triple of smooth functions H a : P — > R one associates a triple of vector fields 
by (no sum on a) X a _l uo a = dH a . The hyperhamiltonian vector field X on P 
associated to the triple {H a } is the sum of these, X := X)a=i X a ', it is trivial to 
check that the X a , and therefore X, are uniquely defined. Each X a is obviously 
Liouville, and so is A. 

On the (p+ 1) dimensional manifold M = RxP (denote by t the coordinate 
on R) the time evolution under A is described by the vector field Z — d t + X. 

Let p a be one-forms (non unique, and possibly defined only locally) satisfying 
dp a = u> a , and ( a the (2N — l)-th exterior power of w a . Define (with s = ±1 
taking care of orientation matters [4] ) 

3 3 
$ = P aA ^ a + ( 6iVs ) ^" Ca A dt . (20) 

QL — 1 QL — 1 

It is immediate to check that dt? is nonsingular, and that Z_ldt9 = 0, 
Z _\dt = 1. It follows from our general discussion that the vector field Z is 
also obtained by a maximal degree variational principles based on the form 1?. 

5 Examples — II 

In this section we consider some concrete and simple - but physically relevant 
- examples: the Euler top, the ABC flow, and the motion of a charged particle 
and of a spin in a magnetic field. 
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5.1 The rigid body 

Let us consider the Euler equations for rotations of a free rigid body around 
its center of mass in three-dimensional space with orthonormal basis vectors 
(e!,e2,e 3 ); we work in the space R 3 of angular velocities, so that x represents 
the angular velocity and x % is the component of the angular velocity in the 
direction e^. The evolution vector field is now given by X — f l di where 

f 1 := mx 2 x 3 ; f 2 := ^ 2 x 3 x 1 ; f 3 := fj, 3 x 1 x 2 

(21) 

Mi = {h - h)/h ; ^2 - (/ 3 - h)/h ; Ms = (h - h)/h • 

The volume form is = da; 1 A dx 2 A dx 3 , hence X _l £1 = f 1 dx 2 A dx 3 - pdx 1 A 
dx 3 + f 3 dx x Ada; 2 . 

Writing the one-form 7 as 7 = A^dx^ , equation (1) reads /' = eijkdjAk or, 
in vector notation, f = rot(A). The solution to this is, up to an exact form 
7' = dg (ja scalar function), 

Ax = \v2x\x 3 ) 2 , A 2 = ^ 3 *V) 2 , A 3 = \^x 3 {x 2 ) 2 . (22) 

Therefore, see (3), the 9 appearing in the variational principle for the Euler rigid 
body equations is (with the A^s given above), 

3 

= x 1 dx 2 A dx 3 - A * ( dxl A d *) • ( 23 ) 

i=l 

5.2 The ABC flow 

The so called ABC flow has been introduced by Henon and is of interest in fluid 
dynamics; it lives in R 3 is characterized by the fact that, with X = f l di, the 
vectors f and rot(f) are collinear (see e.g. [1]). 
We have explicitely 



(24) 



f 1 = Asinix 1 ) + C cos(x 2 ) , f 2 = Bsin^ 1 ) + Acos(x 2 ) 
f 3 = Csin(x 2 ) + 5cos(x 1 ) . 

The volume form is again ft = dx 1 A dx 2 A dx 3 , hence 

X_in = [Csin(x 2 ) + J Bcos(x 1 )]dx 1 A dx 2 

- [Bsin(x 1 ) + Acos(x 3 )]dx 1 A dx 3 (25) 
+ L4sin(x 3 ) + Ccos(x 2 )] dx 2 A dx 3 . 

The solution to X = dj is, up to exact forms 7', given by 

7 = [Asin(x 3 ) + Ccos(x 2 )]dx 1 + [_Bsin(x 1 ) + Acos(x 3 )]dx 2 
+ [Csin(x 2 ) + J Bcos(x 1 )]dx 3 . 

The 1? defining our variational principle for the ABC flow is therefore, up to an 
inessential closed form, given by 

= x 1 dx 2 A dx 3 - 7 A dt . (27) 



(26) 
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5.3 Particle motion in a stationary magnetic field 

We consider now a point particle of mass m and charge q moving in three- 
dimensional space under the effect of magnetic field B (no electric field) . Then 
the equation of motions are x = fc(v x B), with k = q/m and with x the cross 
(vector) product in R 3 . We can rewrite these as a first order system in R 6 as 

v = k (v x B) (28) 

(if B is also constant in space, we can just consider the second equation and 
work in the three-dimensional space of velocities). The corresponding vector 
field, obviously with zero divergence, is thus 

X := v\d/dx l ) + w l {d/dv l ) (29) 

where we have defined, with e the Levi Civita tensor, w % :— ke 1 ^v 3 B l . 

We write u)i :— dx l A dv l (no sum on i), and the volume form £1 in R 6 will 
be written as {I — w\ A u>2 A lo 3 . 

We have, with this notation, 



X = (v^-dv 1 - w'di 1 ) A u>2 A 0J3 + (v 2 dv 2 - w 2 dx 2 ) A 0J3 A + 
+ (v 3 dv 3 - w 3 dx 3 ) Awi A lu 2 ] ■ 

Let us consider a 7 e A 4 (R 6 ) of the form 

7 = 71 - 72 ; 

71 := (1/2) [(v v ) 2 uj2 Auj 3 + (o 2 ) 2 w 3 Aui + (» 3 ) 2 wiAw 2 ] , 

72 := [F a (x)v 2 + F b (x)v 3 ]oj 2 Auj 3 + [G a (x)v 1 + Gb(x)v 3 ] ui 3 A lo\ 

+ [Haix^ 1 + H b (x)v 2 } Wl Aw 2 . 

With these choices, we obviously have that 

d 7l = v l [(d/dx^-lSl] = 

= v 1 (dv 1 A UJ2 A 0J3) + v 2 (dv 2 A UJ3 A UJ\) + v 3 (dv 3 A 0J1 A uj 2 ) 



(30) 



(31) 



(32) 



and thus, recalling that B does not depend on v, X _\VL = d~f if and only if 

(dFa/dx 1 ) = kB 3 (dFb/dx 1 ) = -kB 2 (dG a /dx 2 ) = -kB 3 



(dGb/dx 2 ) = kB 1 (dH a /dx 3 ) = kB 2 {8H b /dx 3 ) = -kB 1 



(33) 



Needless to say, we can always find functions F a , ...,Hb satisfying these, just by 
integrating B l in different variables = 1, 2, 3). 

As for a, we can e.g. choose 

a = (1/3) (x 1 dv 1 A uo 2 Aw 3 + x 2 dv 2 A uj 3 A ui + x 3 dv 3 A wi Aw 2 ) ; (34) 

the form corresponding to X given in (29) is then immediately obtained by 
(3) using this and (31). 
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5.4 Spin motion in a magnetic field. 

As a last example, we consider a particle with spin 1/2 in a spatially constant 
magnetic field, and the evolution of its spin wave function as described by the 
Pauli equation. Disregarding evolution of the wave function associated to the 
spatial coordinates, this reads 



r/vj> 

V = in (B-S)tf • 

dt ' ' 



(35) 



Here n — Airfi/h is a dimensional constant, ^ is a two-components spinor, 
* = (^) + ,^_) with i[)±(i) E C and |*| 2 = 1, the real vector B is the magnetic 
field, with components B(t) = (B x ,B y ,B z ), and S is the vector spin operator 
with components the Pauli a matrices, so that the linear operator M := B • S 
appearing in (35) is given by 



M 



B z B x — iB y 
B x + iB y —B z 



(36) 



The equation (35) can be rewritten as an equation in R rather than in C 2 . 
In order to do so, we rewrite ip± separating their real and imaginary part as 
ip± = \± + K±- With the isomorphism C 1 ~ R 2 given by 1 ~ (1, 0), i ~ (0, 1), 
the operator of multiplication by i is represented in R 2 by the real antisymmetric 
matrix with J12 = — 1 = — J21, and thus 



iM : 



B Z J 

'By I + B X J 



Byl + B X J 

-BJ 



(37) 



Finally, we have obtained that the R 4 representation of equation (31) is 
given by 

f t =*M, (38) 



where 





-B z 


-B x 
-B,, 



By 

B x 


-B z 



By 

B z 
/ 



(39) 



It is easy to check by explicit computation that this equation is hyperhamilto- 
nian (with orientation corresponding to s = — 1), with the symplectic structures 
given by 



uji = dx 1 A dx 3 + dx 2 A dx 4 , oj 2 = dx A A da; 1 + da; 2 A dx 3 
w 3 = dx 2 A da; 1 + da; 3 A dx 4 ; 

and the hamiltonians given by 

HHt,t) = (l/2)B y (t)U\\ 2 , H 2 {t,t) = {l/2)B x {t)U\\ 2 , 
H 3 ^t) = (l/2)B z (t)U\\ 2 ■ 



(40) 



(41) 
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We can choose e.g. 

p 1 = x 1 dx 3 + x 2 dx 4 , p 2 = x 4 dx 1 + x 2 dx 3 , p 3 = x 2 dx 1 + x 3 dx 4 ; (42) 

with these and (40), (41) above, the expression for d is immediately read from 
(20), with s = -1. 
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Appendix. 

Hodge duality and the variational principle 

In this appendix we show how the forms entering in our discussion is related to 
the vector field X via Hodge duality. 

Consider a smooth m-dimensional orientable riemannian manifold (M, go) 
with local coordinates (i = 1, m). In a local chart, vector fields and one- 
forms admit the bases made of 9, := <9/<9£* and d£*; these are dual to each other 
and we write di = d£ 4 and d£* = <3j to denote this duality. We denote the Levi 
Civita (completely antisymmetric) covariant tensor by e Ml ... Mm , with value ±1 
according to the parity of the permutation (/ii, ...fi m ). 

The Hodge star is the linear map * : A r (M) — > A m ~ r (M) defined by 



* (d^ 1 A ... A d^) := V ^ UI ^-^ Vr+1 ... Vm dr r+1 A ... A dC m ■ 

[m r). 

Let, as in the main body of the paper, P be a smooth orientable riemannian 
manifold of dimension A*" with metric g and volume form ft, on which a vector 
field X satisfying (1) for some 7 S A N ~ 2 (P) is defined. We consider the manifold 
M = R x P, with metric g n = S <£> g (so that \g \ = \g\) and coordinates 
(t 7 x , x ). 

Let us now consider the form 1? defining the variational principle for X = 
f l di, see theorem 3 above, and given by (3). We have $ — a + dt A 7 and thus 
di? = ft + dt A d 7 . Recalling (1), we have dtf = ft + f [dt A (d t _l ft)] ; the Hodge 
dual of this is 

N 



*(dtf) = *(ft) + / i *[(0 i _lfi)Ad*] = 



If now we consider the vector field which is dual of this one-form in the sense of 
the duality between forms and vector fields, we have 



*(dt?) - V\g^\ [dt + x] = V\g^\z . 

By reversing the argument, or by acting first with the tilde and then the star 
operators on this equation, we have that: 

Lemma. The form ■& e A N ~ 1 (M) defining the variational principle associated 
to the Liouville vector field X satisfies d$ — y/\g~ 1 \ * (Z). 

Note that this condition completely determines the variational principle: 
indeed it identifies ^ up to a closed form, which has no role in the variation of 
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